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Abstract 

Improved  understanding  of  midwave  infrared  (1-5  micron)  spectral  emissions 
from  detonation  fireballs  is  needed  to  develop  phenomenological  models  for  battle  space 
optical  forensics.  The  ability  to  measure  radiance  over  a  wide  band  pass  at  arbitrary 
resolutions  make  Fourier-transform  spectrometers  (FTS)  an  attractive  tool.  However, 
interferometer  based  spectroscopic  measurements  can  be  corrupted  when  the  observed 
intensity  changes  during  data  acquisition.  While  small,  random  fluctuations  in  scene 
intensity  translate  into  noise,  systematic  variations  introduce  scene-change  artifacts 
(SCAs)  into  Fourier-transformed  spectrum.  Approximating  a  detonation  fireball  as  a 
blackbody  radiator,  modified  with  an  atmospheric  attenuation  ,  with  an  exponentially 

cooling  temperature  T (t)  =  ,  double  sided  interferograms  from  an  ideal  Michelson 

interferometer  were  simulated  and  converted  to  a  measured  spectrum  to  characterize 
the  nature  and  magnitude  of  SCAs.  With  =  1500.^f ,  changing  scene  spectra  on 

0  <  V  <  7900cm"^  at  Sv  =  0.04cm”Vesolution  were  computed  with  percent  changes  in 
T(t)  between  0  and  80%.  The  real  part  of  (Re{L^} )  is  well  approximated  by  the 

instantaneous  spectrum  (L^pjj )  at  zero  path  difference.  Differences  between  Re{L^}  and 
Lzpd  ( -  Lzpp, )  is  largest  when  is  highly  structured  and  these  differences 
increase  with  the  amount  of  temperature  change  during  the  measurement.  Similarly,  the 
imaginary  part  of  (Im{L^} )  exhibits  highly  structured  features  whenever 


rv 


changes  rapidly  with  v  .  In  general,  Im  {4}>  RejL^j-L^p^  agreeing  with  previous 
theoretical  work  and  the  increase  in  magnitude  of  with  temperature  change  is 

more  pronounced  than  Re{L^}  -  .  These  trends  were  independent  of  initial 

temperature,  spectral  resolution,  or  choice  of  .  For  a  MOPD  =  .5cm,  a  4%  reduction 
in  temperature  during  a  single  sean  produced  an  RMS  value  of  5/JW  /  sr  !  cm^  /  cm”'  in 
Im{L^} .  A  21%  temperature  ehange  was  needed  to  produce  the  same  RMS  value  for 
Re  {L^  }  -  L^pjj .  For  comparison,  the  RMS  value  of  was  507 /jW  I  sr  I  cm^  /  cm”'  and 

\25  jJW  I  sr  I  crn^  /  cm”'  for  a  20%  and  60%  temperature  change,  respectively. 

SCAs  were  also  studied  with  an  ABB/Bomem  MR- 154  FTS  by  measuring  a  time- 
varying  broadband  source.  Using  an  Oriel  Q-series  ineandescent  lamp,  the  FTS  operating 
at  Icm  *  resolution  observed  the  changing  spectral  radiance  as  the  lamp’s  input  power 
was  smoothly  varied  during  interferogram  acquisition.  A  2.52%  relative  change  in  ZPD 
voltage  during  the  scan  produced  a  few  spectral  lines  in  Im{L^}  that  were  54.8%  above 
the  instrument  RMS  NESR. 
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SCENE  CHANGE  ARTIEACTS  IN  EOURIER  TRANSEORM  SPECTROSCOPY  OP 
TEMPORALLY  CHANGING  SOURCES 


I.  Introduction 


Background 

Numerous  defense  and  civilian  applications  require  spectral  measurement  of 
temporally  evolving  sources.  Civilian  applications  may  remove  scene  change  effects 
through  repeated  sampling  which,  at  its  simplest,  amounts  to  averaging  out  the  errors. 
Many  defense  applications  can  ill  afford  the  luxury  of  repeated  sampling.  The  events  are 
either  singular  as  is  the  case  with  muzzle  flashes  and  Improvised  Explosive  Devices 
(lEDs),  or  reaction  is  time  critical  as  is  the  case  with  rocket  plumes.  Using  interpreted 
spectra  to  determine  the  nature  of  an  explosive  device  is  a  critical  tool  in  battlespace 
awareness  and  forensics. 

The  Air  Porce  Institute  of  Technology  Remote  Sensing  Group  (APIT  RSG)  is 
pursuing  measurements  of  transient  events  to  develop  phenomenological  models  for 
battlespace  events.  One  tool  in  use  to  acquire  these  signatures  is  Pourier  transform 
spectroscopy  (PTS)  which  allows  the  investigation  of  spectrally  interesting  phenomena. 
The  chemicals  in  explosions  produce  similar  spectra,  but  over  small  frequency  bands, 
unique  features  create  the  possibility  to  investigate  compositions  and  concentrations  to 
better  classify  chemical  reagents.  Most  of  the  bands  of  interest  are  molecular  bands  in 
the  infrared  where  FTS  is  strongly  suited  for  analysis.  FTS  provides  many  benefits  over 
other  spectral  instruments  such  as  radiometers  and  grating  spectrometers.  When  viewing 
transient  events,  the  number  of  photons  reaching  observation  equipment  can  be  small. 
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especially  when  viewing  at  great  distances,  also  making  FTS  more  suitable  than  a  grating 
spectrometer.  FTS  is  able  to  view  large  frequency  ranges  at  high  resolution  where  a 
grating  spectrometer  might  require  multiple  gratings  to  view  all  the  bands  of  interest. 

For  many  systems,  chemical  concentrations  and  temperatures  can  be  revealed 
through  FTS  techniques.  During  many  transient  events  of  interest,  the  scene  changes 
over  a  time  interval  of  data  acquisition  causing  changes  in  interferogram  intensity  not 
caused  by  an  optical  path  length  difference.  The  nature  of  quickly  changing  scenes  on 
FTS  calculated  spectra  is  not  well  understood  and  needs  to  be  characterized  for  correct 
scene  interpretation  to  exploit  transient  event  signatures.  Using  FTS  to  acquire  transient 
event  signatures  might  seem  counterintuitive  because  of  their  rapid  nature;  however, 
previous  observations  from  AFIT  RSG  of  transient  events  have  captured  meaningful 
spectra^  The  extent  to  which  the  spectra  are  corrupted  needs  to  be  investigated  to 
determine  when  certain  approximations  and  certain  conditions  can  provide  accurate  FTS 
spectra  even  when  the  scene  is  changing  rapidly. 

Objectives 

This  research  intends  to  describe  the  effect  of  exponential  temperature  decays  in 
black  body  radiance  profiles  on  spectra  produced  through  an  ideal  FTS  system  and  to 
determine  the  validity  of  using  the  zero  path  difference  (ZPD)  spectrum,  the 
instantaneous  spectrum  corresponding  to  when  the  Michelson  is  at  ZPD,  as  an 
approximation  of  the  calculated  real  spectrum  for  an  exponentially  cooling  black  body. 
Also  whether  the  calculated  real  spectrum  corruption  can  be  inferred  based  upon  the 
imaginary  spectrum  and  how  it  relates  to  the  calculated  real  spectrum  as  a  function  of 
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temperature  change.  This  will  help  determine  the  conditions  when  FTS  is  corrupted 
through  scene  change  and  when  it  is  viable  for  observing  such  transient  events.  This 
work  primarily  looks  at  RMS  values  to  quantitatively  look  at  entire  spectra.  A  more  in 
depth  look  of  specific  spectral  features  would  also  be  beneficial  for  analysis  of  spectral 
features  used  in  forensic  analysis. 
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II.  Theory 


FTS  background  and  introduction  to  the  complex  spectrum 

A  Michelson  interferometer  splits  incident  light  into  two  paths.  Within  one  arm 
of  the  interferometer,  the  path  length  is  changed  and  when  the  light  is  recombined, 
interference  occurs.  The  intensity  of  the  recombined  light  is  measured  as  a  function  of 
optical  path  difference  (OPD)  and  the  intensity  of  the  interferogram  given  by  Brault  will 
follow'^^^ 

^00  .  ^00 

=  jg  ri{v)ri,niy)~  L{v)dv  +  ^^  ri{v)L{v)co^{ln:xv)dv ,  (1) 

where 

£  V(y)Vm(yy^L{v)dv  (2) 

is  the  DC  component  of  the  interferogram  and 

poo 

ri{v)L{v)cos{2Kxv)dv  (3) 

Jo 

is  the  AC  component. 

rjiy)  is  the  interferometer  efficiency  and  is  given  by  the  expression 

(4) 

rj^iv)  is  the  beam  splitter  efficiency  and  is  given  by  the  expression 

%{v)  =  AR{v)T{y)  (5) 

where  R{v)  is  the  wavenumber  dependent  reflectance  of  the  beam  splitter  and  T{v)  is  the 
wavenumber  dependent  transmittance  of  the  beam  splitter.  For  an  ideal  case,  half  of  the 
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incident  light  at  every  wavenumber  is  reflected  and  half  is  transmitted  resulting  in 

77^  (v)  =  ^  .  (i7)  is  the  efficiency  of  all  the  other  optics  and  would  be  unity  in  the  ideal 

case.  77„(v)  is  the  modulator  efficiency  which  is  the  ZPD  amplitude  divided  by  the 
averaged  baseline  modulation  amplitude  of  the  interferogram.  The  ideal  case  also 
ignores  the  wavenumber  dependent  phase  shift,  (f>- ,  eaused  by  all  opties  in  the 
interferometer.  The  AC  component  of  the  interferogram  ean  be  expressed  as 

poo 

/(x)  =  ri{v)L{v)co?,{27rvx-^-)dv .  (6) 

J  0 

The  remaining  components  creating  the  interferogram  are  L(v) ,  the  incident  spectral 
radiance,  and  x ,  the  OPD.  Only  the  AC  component  is  needed  at  this  point  to  continue. 

From  simplifications  found  in  Sromovsky,^^^  a  complex  spectrum  G(i7)  defined  for 
all  V  can  be  introduced  as 

^77(v)L(F)exp(/^2i-),v>0 

G(v)=  ^  .  (7) 

-  rii-^)L{-v)  exp(-/^?i_- ),  v  <  0 

The  real  transform  of  the  complex  spectrum  will  yield  the  interferogram  shown  by 

poo 

/(x)=  G{v)Qx^{-27rivx)dv .  (8) 

J-00 

Consequently,  the  speetrum  is  now  given  by  the  eomplex  Fourier  transform  of  the  real 
interferogram'^'^l  The  original  spectrum  can  now  be  defined  in  terms  of  this  transform  as 

2  (“» 

L(i7)  = - exp(-i^i )  7(x)  exp(2;Tixv)dx .  (9) 

7(v) 
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For  a  typical  Michelson  interferometer  which  produces  interferograms  that  are 
discretely  sampled  at  uniform  OPDs,  a  discrete  Fourier  transform  best  represents  the 
spectrum.  The  limits  of  integration  with  respect  to  OPD  are  finite  and  the  Fourier 
transform  defining  the  spectrum  is  best  approximated  through  the  finite  series 

"2,  ^  A^/2-l 

^  /(xjexp(2;r/vxj  (10) 

k=-N/2 

where  X  is  maximum  OPD  (MOPD)  and  N  is  the  number  of  points  in  the  interferogram. 

I  I  N 

If  |v|  <  ,  then  the  spectrum  is  sampled  adequately  and  can  be  represented  by  only  N 

wavenumbers  given  by 


2Y  N/2-1 

H(vJ  =  —  ^  I(Xf^)exp(27rink/  N) 

^  k=-N/2 


(11) 


where 


v„  =■ 


2X 


(12) 


and 


n  =  -A^/2,-A^/2+l,...,A^/2-l.  (13) 

The  preceding  development  assumes  a  static  source,  L(v) ;  one  that  is  constant 
during  the  measurement.  If  the  radiance  of  the  incident  light  changes  while  the  device  is 
scanning,  scene  change  artifacts  (SCAs)  in  the  spectrum  result  from  changes  in 
interferogram  not  caused  by  a  change  in  optical  path  length  difference. 
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Previous  theoretical  scene  change  work 

The  nature  of  SCAs  can  be  better  understood  by  considering  smooth  time  varying 
sources.  Following  previous  work  from  Kick  et  al.^^^  and  using  an  ideal  Michelson 
interferometer  so  there  are  no  phase  errors,  if  a  scene  were  changing  linearly  with  time, 
shown  in  Figure  1,  the  spectrum  could  be  given  by 

L{t,v)  =  mv)  +—{L,iy)-my)) 

hot 

where  is  the  total  measurement  time,  Lj  (i7)  is  the  spectral  radiance  at  the  beginning  of 
the  measurement,  and  is  the  spectral  radiance  at  the  end  of  the  measurement.  As 
time  passes,  the  radiance  at  each  wavenumber  changes  linearly  fromL|(i7)  toL3(i7)  until 


(14) 
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Figure  1.  Linear  (blue)  and  quadratic  (green)  variations  in  radiance  with  time  for  a 

single  wavenumber. 

With  a  change  in  variables  from  time  to  path  difference,  the  spectral  radiance  as  a 
function  of  wavenumber  and  path  difference  is  given  by 

L(x,  V)  =  ^  (p)  _  (P))  (15) 

2  Ax 

where  Ax  is  the  total  distance  the  mirror  scans  and  x  is  the  instantaneous  path  difference 
Ajc  At 

and  ranges  from  to  —  .  The  interferogram  from  a  continually  sampling 

spectrometer  would  then  be  the  inverse  Fourier  transform  of  Equation  (15)  given  as 

/(x)  =  3-'{L(x,v)}  =  ^^1^^^^^  +  ^(/3(x)-/,(x)).  (16) 

2  Ax^  ^ 
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This  is  the  interferogram  captured  by  the  Michelson.  To  calculate  the  spectrum  the 
Fourier  transform  of  the  interferogram  is  needed.  Kick  et  al.  show  that  if  Equation  (16)  is 
transformed  back  into  a  wavenumber  domain,  the  resulting  effect  on  the  spectrum  is 
shown  by^^^ 

L  (vi  A(0+4(0  ,  i  a(4<0-A<0)  (17. 

"  2  Intix  dv  ' 

The  imaginary  part  comes  from  the  Fourier  transform 

3{^(/3(x)-/j(x))}  (18) 

which  uses  the  Fourier  transform  relationship 
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{x"/W}  = 


A" 


d-fiv) 

dv" 


Knowing  a  linear  scene  change  provides  an  imaginary  theoretical  component 

i  d(L,(v)-L,(v)) 


Yi  =  lm(L(v))  = 


IttAx 


dv 


the  imaginary  part  of  the  spectrum  can  be  used  to  distinguished  SCAs. 

Similarly,  the  same  analysis  can  be  used  for  a  scene  radiance  that  is  changing 

quadratically  with  time,  shown  in  Figure  1,  resulting  in 

f  m  ’  a(2i,<r)-42^„(r)  +  2Z,(P))  i  a(i,(r-)-Z,(.:)) 

'  4^'- Ax'-  dv^  iTdxx  dv 

where  L^pj^  corresponds  to  the  ZPD  spectrum  observed  by  the  spectrometer^^l  There  now 

is  a  second  order  correction  to  the  observed  real  spectrum  given  by 

1  d{lL,{y)-AL^p^{v)  +  lU,{v)) 


Yr  =  — 


dv^ 


(19) 


(20) 


(21) 


(22) 
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and  the  first  order  correction  to  the  imaginary  piece  remains  the  same  as  the  linear  case. 

The  same  analysis  can  be  used  with  higher  order  terms  for  more  accurate  approximations. 

Connection  between  previous  work  and  discrete  FFT 

This  analysis  is  not  obviously  applicable  to  the  discrete  Fourier  transform 
described  above  and  therefore  a  real  Michelson  system.  In  the  discrete  transform,  the 
OPD  does  not  range  from  ^aDtooo,  which  is  necessary  for  the  full  transform  found  in 
Kick  et  al.  To  limit  the  extent  of  the  integral  and  to  apply  their  analysis  to  discrete 
transforms,  let  the  interferogram  intensity  equal  the  convolution  given  by 

i(x)  =  I(x)  *  rect{x,  MOPD) .  (23) 

This  will  result  in  the  Fourier  transform  giving  the  spectrum  as 

3  =  j"  l{x)dx  (24). 

However,  Im(x)  =  0  when  |x|  >  MOPD  resulting  in 

The  spectrum  is  not  continuously  sampled,  therefore  the  spectrum  given  by  Equation  (25) 
is  discretized  and  results  in  Equation  (10).  This  allows  for  the  discrete  sampling  and 
finite  MOPD  characteristics  of  a  Michelson  to  be  represented. 

Connection  to  detonation  fireball 

Within  the  realm  for  detonation  fireballs,  exponentially  decaying  models  can  accurately 
describe  the  radiance  change  over  time.  Defining  a  decay  constant  that  controls  what 
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temperature  change  occurs  creating  a  smoothly  varying  source  over  a  half  scan  period  for 


a  Michelson  allows  for  comparisons  to  Kick  et  al.  analysis.  When  the  temperature 
difference  is  small  enough  so  the  exponential  decay  creates  a  spectrum  that  can 
accurately  be  described  as  a  linear  or  quadratic  decay,  the  real  and  imaginary  spectrum 
can  accurately  be  described  by  the  mathematics  presented  in  Kick  et  al.  Regardless  of  the 
mathematical  representation  of  SCAs,  the  magnitude  of  these  features  can  be  compared 
to  a  noise  equivalent  spectral  radiance,  NESR,  to  see  if  scene  change  creates  artifacts  that 
rise  above  the  instrument  noise  floor.  This  will  determine  the  corruptness  of  using  FTS 
and  if  no  SCAs  are  above  the  inherent  noise  floor  of  the  system,  there  is  no  corruption 
from  SCAs. 
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III.  Modeling 


Model  introduction 

An  ideal  detonation  fireball  is  modeled  as  a  black  body  radiator  with  an 
exponentially  decaying  temperature  given  by 

B{v,T{t))  =  Ihc^v^  -  (26) 

ehn,) 

where  temperature  follows 

T{t)  =  T^e-'^‘  (27) 

and  Tq  is  the  initial  temperature  and  ^  is  a  decay  constant.  Including  atmospheric 
transmission,  r^,„(v) ,  resulted  in  a  time-evolving  expression  for  radiance  as 

=  (28) 

By  inputting  a  range  of  optical  path  length  differences  based  on  a  HeNe 
wavelength  sampling,  632.816  nm,  an  interferogram  was  created  by  using  the  cosine 
transform 

poo 

y(v.)  =  D{v)r]{y)L{t^,v)cos{2mx^v-<l)-)dv  (29) 

where  D{v)  is  the  detector  response  at  a  given  wavenumber  and  for  the  model  is  unity 
over  the  wavenumbers  from  1800  cm  *  to  7900  cm  *  and  zero  elsewhere.  is  the 
wavenumber  dependent  phase  shift  and  in  the  model  was  the  ideal  case  of  zero.  x.  is  the 

optical  path  length  difference  of  the  two  interferometer  arms  corresponding  to  the  time 
interval,  v  is  the  wavenumbers  over  which  the  model  is  made,  and  L(t.,v)  is  the 
spectrum  created  above  at  a  specific  time  matching  the  measurement  being  made  while 
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OPD  =  X. .  At  each  successive  HcNe  wavelength,  the  spectram  has  evolved  one  time 
step,  so  unique  values  of  L(t^,v)  corresponds  to  each  OPD  value,  x. ,  through  the  relation 

X.  =  vt.  (30) 

where  v  is  the  mirror  velocity. 

The  model  was  created  over  the  wavenumber  range  from  0  cm'^  to  7901  cm'\  the 
Nyquist  frequency  corresponding  to  the  HeNe  wavelength  sampling.  By  allowing  the 
temperature  to  decay  to  20  percent  of  its  initial  temperature,  the  decay  constant  k  could 
be  solved  for  by  knowing  the  number  of  points  in  the  interferogram  which  is  equal  to  the 
total  number  of  time  steps  taken. 

The  number  of  samples  of  the  frequency  axis  for  a  given  MOPD  must  insure  the 

I  I  N 

condition  |v|  <  is  satisfied.  If  the  frequency  axis  is  not  adequately  sampled,  artifacts 

not  caused  by  scene  change  can  contribute  to  the  spectrum.  Solving  for  the  needed 
spacing  in  the  frequency  axis  gives  the  relationship  for  the  number  of  points  as 

=  Nyquist{2MOPD)  (3 1 ) 

but  more  points  are  needed  to  define  the  derivatives  at  each  spectral  coordinate,  so  using 
25  times  the  output  of  Equation  (31)  was  more  than  sufficient  to  accurately  create  the 
derivatives  needed  to  mathematically  represent  scene  change  in  Equation  (21). 

The  interferogram  created  through  this  model  will  look  like  Figure  2  and 
depending  on  how  much  the  black  body  is  changing  in  temperature,  the  baseline 
modulation  amplitude  at  the  beginning,  on  the  left,  will  be  greater  than  that  of  the  end,  on 
the  right,  because  the  intensity  of  the  spectrum  is  greater  at  that  time. 
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Figure  2.  Interferogram  of  a  black  body  cooling  during  data  acquisition.  MOPD= 
.25  cm.  Initial  Temperature  =  1500  K.  Final  Temperature  =  600  K. 

While  there  is  a  DC  component  for  the  transform  because  the  spectrum  is 
evolving,  the  AC  coupling  of  the  Michelson  system  is  such  that  the  DC  component  is 
subtracted  faster  than  the  modulated  part  of  the  interferogram  resulting  in  an 
interferogram  with  a  baseline  of  zero. 

The  interferogram  is  transformed  through  a  discrete  Fourier  transform,  given  by 
Equation  (1 1),  to  give  the  resulting  spectrum  and  the  real  and  imaginary  components  are 
compared  to  the  pieces  found  in  the  Kick  et  al.  analysis.  The  spectrum  corresponding  to 
the  interferogram  in  Figure  2  is  presented  in  Figure  3.  The  spectral  features  are  the 
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atmospheric  transmission  features  over  the  given  wavenumber  region  for  an  atmospheric 
path  length  of  0.5  km.  The  process  of  turning  the  interferogram  into  a  spectrum  also 
cancels  any  effects  given  by  the  previous  processes.  From  Equation  (9),  a  factor  of  a  half 
from  the  system  efficiency  factor,  rjiy)  ,  and  another  factor  of  two  from  the  conservation 
of  energy  in  defining  a  symmetric  spectrum  defined  over  all  v  in  Equation  (7)  results  in  a 
factor  of  four  being  applied  to  the  integral  defining  the  spectrum  after  the  Fourier 
transform  to  undo  the  effects  of  the  interferometer. 


Figure  3.  Calculated  spectrum  from  interferogram  presented  in  Figure  2. 
Imaginary  spectrum  is  on  top  and  real  spectrum  is  on  bottom. 
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The  imaginary  structures  correspond  to  the  regions  in  the  real  spectrum  where  there  are 
abundant  spectral  features  relating  to  a  highly  structured  . 

Using  this  model  will  allow  the  simulation  of  a  Michelson  observing  an 
exponentially  cooling  black  body  through  atmosphere  over  0  to  7901  cm'\  The 
calculated  spectrum  can  be  compared  to  the  Kick  et  al.  analysis.  By  separating  the  real 
and  imaginary  components  and  calculating  the  theoretical  components  through  the 
derivative  found  in  Equation  (21),  the  implementation  of  Kick  et  al.  can  be  verified. 

It  was  unknown  as  to  whether  the  mathematical  treatment  in  Kick  et  al.  accurately 
accounts  for  the  spectral  information  contained  in  the  time  evolving  nature  of  scene 
change.  After  taking  a  Fourier  transform,  this  would  create  spectral  information  not 
attributed  to  the  beginning  spectrum.  To  address  the  problem  of  aliasing,  oversampling 

HeNe 

was  used.  With  a  wavelength  as  small  as  126.4  nm,  — ^ — ,  no  spectral  features  were 

observed  outside  of  the  black  body  curve  up  to  the  Nyquist  frequency  and  no  noticeable 
changes  to  the  spectrum  occurred  indicating  additional  spectral  information  was  not 
present. 

Metrics 

By  modeling  a  rapidly  cooling  Planckian,  the  nature  of  SCAs  can  now  be 
examined.  Specifically,  the  measured  spectrum  L^(v')  will  be  compared  to  the 

instantaneous  spectrum  at  ZPD.  For  several  comparisons,  root-mean  squared,  RMS, 
values  were  used  to  compare  trends  and  entire  spectra. 
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SCA  metrics  need  to  be  created  to  help  understand  the  extent  to  which  spectral 
measurements  of  evolving  sources  become  corrupted  when  using  FTS.  Moreover,  being 
able  to  determine  if  scene  change  is  significantly  corrupting  the  real  spectrum,  Re{Lj^}  , 

by  examining  the  imaginary  spectrum,  Im{L^} , would  be  very  beneficial.  A  useful 
relationship  is  comparing  the  imaginary  spectrum  to  the  corruption  in  the  real  spectrum 
given  by  the  difference  between  Re{Lj^}  and  ,  Re{L^}  -  L^pj^ .  Showing  the  change 

in  the  RMS  value  of  Im{L^}  along  with  the  RMS  value  of  Re{L^}  as  both  evolve 

as  a  function  of  temperature  change  shows  how  the  connection  between  increased 
imaginary  components  and  the  deviation  from  the  L^pp, .  More  specifically,  the  real  SCA 
RMS  magnitude  given  by 

SCA^^=RMS(R^{Lj-L,pp,)/RMS(Re{Lj)  (32) 

and  the  imaginary  SCA  RMS  magnitude 

=RM5(Im{L„})/RMS(Re{Lj)  (33) 

are  used  to  describe  the  magnitude  of  SCAs  as  they  change  with  certain  parameters. 

Similarly,  the  maximum  real  SCA  magnitude 

MSC\^  =  Max(Re{L^ }  -  4^^ )  /  RMS  (Re  {4 })  (34) 

and  the  maximum  imaginary  SCA  magnitude 

MSCA^  =Max(Im{4})/RMS(Re{4})  (35) 

are  used  for  analysis. 
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For  different  decay  rates,  the  exponential  curve  creates  a  change  in  the  spectrum 
that  can  be  approximated  by  a  first  or  second  degree  polynomial.  How  well  this 
approximation  fits  the  actual  spectrum  is  seen  in  the  deviation  of  Re{Lj^}  fromL^^^ . 

Comparing  the  difference  between  the  theoretical  predictions  based  on  Kick  et  al.  and  the 
actual  calculated  real  spectrum  using  an  exponentially  cooling  blackbody  will  ultimately 
show  how  well  and  for  what  conditions  the  theoretical  predictions  are  a  good 
mathematical  approximation.  The  trends  of  the  model  will  be  spectrum  dependent, 
however,  by  comparing  different  temperature  changes  and  different  atmospheric 
transmission  features,  a  family  of  curves  can  be  seen  with  the  same  general  trend.  It  is 
also  important  to  consider  the  magnitude  of  SCAs  and  how  they  compare  to  a  detection 
instrument’s  noise  equivalent  spectral  radiance,  NESR. 

Model  analysis 

For  Figure  4,  the  exponential  temperature  decay  from  1500  K  to  1200  K  creates  a 
change  in  the  spectrum  that  can  be  closely  approximated  by  a  quadratic  polynomial  at 
every  wavenumber  resulting  in  nearly  all  lines  to  be  indistinguishable.  This  allows  for 
the  mathematical  representation  found  in  Equation  (21)  to  accurately  represent  scene 
change.  Re  {L^  }  closely  resembles  with  the  difference  being  the  second  order 

correction  found  in  the  quadratic  case,  Yr. 
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wavenumber  (cm‘^) 


Figure  4.  FTS  spectrum  of  black  body  with  atmospheric  transmission  path  length  of 
.5  km  and  a  temperature  change  from  1500  K  to  1200  K  during  data  acquisition. 

MOPD  =  .25cm. 

The  corraption  Re  {L^  }  -  is  a  maximum  of  20 /lM  I  sr  I  crn^  /  cm”'  whereas 

Re{L^}  has  a  maximum  near  1400//VT  /  sr  !  cm^  /  cm"'  meaning  the  corruption  is  a  small 

fraction  of  the  total  signal.  The  matching  of  SCAs  to  wavenumber  regions  in  the 
spectrum  near  highly  varying  spectral  features  confirms  the  derivative  relationship  to  the 
magnitude  of  SCAs  found  in  Equation  (21).  Large  values  of  the  derivative  of  the 
spectrum  match  with  increased  magnitude  in  SCAs. 
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Figure  5.  Top)  Difference  between  Yi  and  Im{L^}  for  Figure  4.  Bottom)  Difference 
between  Re {L^ }  -  and  the  Yr  for  Figure  4. 

Figure  5  shows  the  deviation  of  the  calculated  spectrum  from  the  theoretical 
spectrum  for  Figure  4  and  ultimately  how  accurately  Equation  (21)  can  represent  the 
scene  change  effects.  There  is  a  numerical  noise  floor  that  the  features  must  rise  above 
seen  in  0  to  2000  cm'^  on  the  bottom  graph.  For  the  temperature  change  from  1500  K  to 
1200  K  the  mathematical  expressions  accurately  approximate  SC  As  to  within  a  1  to  3 
fjW  I  sr  I  cm^  /  cnf^  in  the  real  spectrum  and  10  to  30 /jW  /  sr  /  crn^  /  cm~^  in  the 
imaginary  spectrum. 
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wavenumber  (cm"') 


Figure  6.  FTS  spectrum  of  black  body  with  atmospheric  transmission  path  length  of 
.5  km  and  a  temperature  change  from  1500  K  to  600  K  during  data  acquisition. 

MOPD  =  .25cm 

For  the  temperature  change  1500  K  to  600  K  found  in  Figure  6,  the  same  analysis 
is  no  longer  accurate  in  mathematically  describing  the  spectral  change  over  time. 

Re{L^}  now  has  a  maximum  value  near  100//VF  !  sr  !  cm^  /  cm“^  which  is  of  equal 

magnitude  to  the  maximum  values  of  Re{L^}  at  450//1F  /  sr  /  cm^  I  cm~^  meaning  the 

corruption  represents  a  larger  fraction  of  the  total  spectral  radiance.  This  is  also  seen  the 
difference  between  maximum  corruption  magnitudes  of  approximately 
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100//VK  !  sr  !  cm^  /  cm  ^  for  Figure  6  and  approximately  20 juW  /  sr  /  cm^  /  cm  ^  for  Figure 
4. 


Figure?.  Top)  Difference  between  Yi  and  Im{L^}for  Figure  6.  Bottom)  Difference 
between  Re  {L^  }  -  and  tbe  Yr  for  Figure  6. 

There  is  significant  difference  between  Re  {L^ }  -  L^p^,  that  cannot  completely  be 
accounted  for  by  Yr.  Similarly,  Im{L^}  begins  to  deviate  from  the  theoretical  first  order 
prediction,  Yi.  Higher  order  terms  must  be  applied  to  account  for  Re{L^}  -  L^pp,  and  the 
difference  between  Im{L^}  and  Yi.  For  the  temperature  change  from  1500  K  to  600  K, 
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the  mathematical  expression  accurately  approximate  SCAs  in  Figure  7  with  a  maximum 
difference  of  40//VF  !  sr!  cm^  /  cw“'  in  the  real  spectrum  and  200//1F  !  sr!  cm^  /  cm“^  in 
the  imaginary  spectrum.  The  magnitudes  of  the  features  are  greater  than  that  of  Figure  5 
due  to  the  greater  temperature  change. 

To  show  how  much  deviation  between  Re  {L^}  and  exists  for  the  total 

spectrum  and  how  it  relates  to  Im{L^  } ,  Figure  8  compares  the  change  in  the  RMS  of 

Im{L^}  with  the  RMS  value  of  Re{L^}  -L^^^both  as  a  function  of  the  exponential  decay 

constant,  k  ,  multiplied  by  the  measurement  time,  .  All  values  are  normalized  by  the 

RMS  of  Re{L^} .  Some  differences  between  MOPD  models  could  be  created  from 

defining  the  number  of  points  in  Equation  (31)  which  depends  on  the  MOPD  and 
ultimately  affects  the  spectral  resolution  of  the  input  spectrum. 
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Figure  8.  SCArc  (solid)  and  SCAim  (dashed)  relationship  with  k*ttot .  Atmospheric 

transmission  path  length  =  .5  km. 

To  better  allow  for  comparisons,  the  x-axis  in  Figure  9  was  transformed  into  a 
percent  change  in  temperature  to  better  show  how  the  RMS  of  Im{L^}  and  the  RMS  of 

Re  {L^  }  changed.  This  also  allows  for  further  comparison  between  different 
spectra.  For  these  spectra,  temperature  changes  of  55%  yield  a  value  of  under 
10%,  meaning  that  for  temperature  changes  of  up  to  55%  the  RMS  of  the  corruption 
remains  under  10%  of  the  RMS  spectrum  value.  If  Re  {L^}  matched  ,  then  the  RMS 

value  of  the  difference  would  be  zero.  In  general,  Im  {4}>Re{4}-L2p^  for  all 
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temperature  changes  and  Im{L^}  grew  faster  than  Re{L^}  with  temperature 
change. 


Figure  9.  SCArc  (solid)  and  SCAim  (dashed)  relationship  with  temperature  change. 

Atmospheric  transmission  path  length  =  .5  km. 

Instead  of  comparing  RMS  values,  the  maximum  values  were  explored  in  Figure 
10.  Sharp  peaks  in  the  imaginary  spectrum  come  from  large  values  in  the  derivative  of 
the  spectrum  which  are  often  the  result  of  spectral  absorption  features.  How  well  the 
theoretical  predictions  match  the  spectrum  at  those  points  is  very  important  to  accurately 
interpret  the  forensics  of  an  explosion.  Certain  spectral  features  can  often  be  several 
times  the  RMS  value  of  the  spectrum. 
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Figure  10.  MSCArc  (solid)  and  MSCAim  (dashed)  relationship  with  k*ttot  on  the  top 
graph  and  with  temperature  change  on  the  bottom  graph.  Atmospheric 
transmission  path  length  =  .5  km. 

While  this  analysis  is  spectram  dependent.  Figure  1 1  and  Figure  12  show  that  by 
changing  the  temperature  range  and  changing  the  atmospheric  transmission  features  as  a 
function  of  path  length  does  not  significantly  change  the  deviation  from  4pd  • 

Figure  1 1  provides  three  different  maximum  temperatures  for  starting  decay 
parameters.  For  the  three  cases  of  starting  temperatures  4000  K,  2500  K,  and  1500  K 
decreasing  to  20%  of  the  original  value,  the  general  trend  is  the  same  showing  that  at 
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percent  temperature  changes  under  50%,  the  RMS  of  Re{L^}  remains  under  10% 

of  the  RMS  value  ofRejL^} .  The  SC  A  RMS  magnitudes  of  the  different  temperatures 

are  in  part  due  to  their  black  body’s  overlap  with  the  wavenumber  region  of  interest  0  to 
7901cm'^  and  how  much  of  the  black  body’s  power  is  contained  within  that  region. 


Percent  Temperature  Change 


Figure  11.  The  three  scenarios  correspond  to  different  initial  temperatures  of  1500 
k,  3000  k,  and  4500  K.  Dashed  lines  represent  the  SCAim  while  the  solid  lines 
represent  the  SCArb.  MOPD  =  .5  cm  and  Atmospheric  transmission  path  length  = 

.5  km. 
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Figure  12  shows  the  same  general  conditions  of  initial  temperature  of  1500  K 
and  an  MOPD  equaling  .5  cm,  only  this  time  with  different  atmospheric  path  lengths  of  .5 
km,  1  km  and  2  km.  This  results  in  the  transmission  function  having  different 
magnitudes  on  certain  spectral  features.  For  all  the  cases,  the  general  trend  holds  and  for 
percent  temperature  changes  under  55%,  the  RMS  of  Re{L^}  - remains  under  10% 

of  the  RMS  value  of  Re{L^}  and  Im{L,,j}  >  Re{L,^}  . 
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Figure  12.  The  three  scenarios  correspond  to  different  atmospheric  transmission 
profiles  as  a  function  of  the  path  lengths  of  .5  km,  1  km,  and  2  km.  Dashed  lines 
represent  the  SCAim  while  the  solid  lines  represent  the  SCArb.  MOPD  =  .5  cm. 
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Figure  13.  Relationship  between  SCAim  and  SCArc.  Axes  are  shown  as  a  percent  of 

the  RMS  of  Re  {4}. 


To  be  able  to  use  the  imaginary  spectrum  to  help  predict  discrepancies  between 
Re {4}  and  4pd  »  Figure  13  shows  their  relationship.  After  the  RMS  of  Im {4}  reaches 

15%  of  the  RMS  ofRe{4} ,  the  RMS  of  Re{4}  ~4rahas  a  linear  relationship  with  the 
RMS  ofIm{4}  •  It  appears  that  SCAjm  is  about  10%  higher  than  SCAre,  meaning  a 
spectrum  with  SCAj^  =  40%  of  Re {4}  has  a  SC\^  =  30%  of  Re {4}  • 
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The  usefulness  of  Equation  (21)  as  an  approximation  for  a  exponentially  cooling 
Planckian  radiance  viewed  by  a  FTS  can  be  seen  in  Figure  14  that  show  the  deviation  of 
the  calculated  values  in  the  model  from  the  theoretical  values. 


Percent  Temperature  Change 


Figure  14.  Dashed  lines  represent  theoretical  predictions  based  on  implementing 
Equation  (21)  based  on  Kick  et  al.  The  solid  line  represents  the  actual  modeled 


values.  Top)  Deviation  of  RMS  values  of  Im{L^}  and  Yi.  Middle)  Deviation  of  the 
RMS  of  Re{L^}  -L^p^and  Yr.  Bottom)  Deviation  of  the  RMS  ofRe{L^}from  the 


RMS  ofL^p^ . 
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The  deviation  from  theoretical  predictions  is  clearly  seen  in  the  graphs  on  Figure 
14.  For  temperature  changes  up  to  40%,  the  values  of  Equation  (21)  provide  an  accurate 
approximation.  For  the  bottom  graph  representing  the  difference  between  the  RMS  of 
L^pjj  and  the  RMS  of  Re{L^} ,  all  lines  lie  on  top  of  one  another  indicating  that  the 


average  values  do  no  significantly  diverge  from  theoretical  calculations  for  RMS  values. 


Percent  Temperature  Change 

Figure  15.  RMS  of  Im{L^}  (dashed)  and  the  RMS  of  Re{L^}  (solid)  values  as 
they  change  with  temperature  change.  The  black  horizontal  line  represents  a  NESR 
ofS/M!  srlcrn^  Ian  ^  that  both  Im{L^}  and  RejL^I-L^p^  must  rise  above  for 

SCAs  to  be  distinguishable. 
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Figure  15  shows  the  evolution  of  RMS  SC  A  magnitudes  with  temperature  change 


and  provides  a  comparison  to  a  NESR  that  a  SCA  magnitude  must  surpass  to  rise  above 
the  noise  floor.  While  the  RMS  values  may  be  under  or  over  the  NESR  individual  SCA 
magnitudes  at  a  specific  wavenumber  may  or  may  not  be  above  the  NESR  indicating  a 
need  for  analysis  of  individual  wavenumber  SC  As. 
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IV.  Experiment 


Experimental  setup 

To  experimentally  characterize  spectral  artifacts  due  to  a  changing  source,  a 
variable  input  lamp  was  observed  by  an  FTIR.  The  ABB-Bomem  MR154,  a  dual  channel 
interferometer  with  an  internal  cold  reference,  was  used  to  collect  spectra  in  the  mid 
infrared.  The  instrument  was  configured  with  InSb  (1.6-5. 6pm)  and  HgCdTe  (2.0pm  - 
20pm)  detectors.  All  scans  with  the  FTIR  were  done  at  a  resolution  of  1  cm  *  to  slow 
down  the  acquisition  rate  resulting  in  running  a  collection  frequency  of  1.17  Hz.  An 
Oriel  Q  series  incandescent  lamp  was  used  as  the  variable  light  source  despite  the  limited 
overlap  in  the  infrared  region  of  the  detectors.  The  lamp  envelope  and  optics  have  a 
passband  of  approximately  0.2pm  to  2.5pm.  Convolving  the  passband  of  the  lamp 
envelope,  the  Germanium  blocking  filter  and  the  InSb  detector  response  yields  an 
approximate  range  of  1.8pm-2.8  pm  for  the  response  due  to  the  lamp  filament.  Only  the 
InSb  channel  was  able  to  capture  the  signal  as  the  HgCdTe  was  outside  the  passband. 

The  setup  in  Figure  16  shows  the  relationship  between  the  lamp  and  the  FTIR. 

All  incident  light  goes  through  the  Germanium  window  to  control  the  wavelengths 
involved  in  the  interferogram  and  to  stop  aliasing  from  higher  wavelengths.  Possible 
scattered  light  can  reach  the  detectors  from  dust  on  the  optics  although  this  is  assumed  to 
be  negligible. 


33 


Collimated  lamp  with 
quart:  lens 


Germanium 

Window 


Bomem  FTIR 


Figure  16.  Experimental  Setup 

For  reference,  30  interferograms  were  taken  for  static  lamp  powers  descending 
from  210  W  in  7  Watt  increments.  A  total  of  ten  measurements  were  made  for  each 
point.  Each  measurement  made  was  a  1/2  scan  resulting  in  each  power  level  being 
scanned  in  two  different  directions.  Each  scan  direction  might  not  have  the  same 
properties  so  they  must  be  accounted  for  separately.  This  gave  five  measurements  in 
each  scan  direction  that  was  averaged  together  to  increase  signal  to  noise. 

Changing  source  interferograms  were  taken  during  power  decay  over  the  same 
Watt  range  as  above;  however,  interferograms  were  taken  while  the  source  was  changing. 
Two  trials  were  performed,  the  first  with  a  total  of  50  measurements,  25  in  each  direction, 
and  the  second  with  75  measurements,  38  in  direction  one  and  37  in  direction  two.  The 
FTIR  started  taking  measurements  and  the  lamp  was  then  triggered  after  some  time  to 
decrease  from  210  W  to  zero.  The  interferograms  were  the  converted  into  spectra  and 
phase  corrected  using  code  provided  by  Bomem. 


Experiment  analysis 

Viewing  the  changing  lamp  provided  insight  into  the  detection  of  SCAs.  The 
imaginary  part  created  from  SCAs  must  be  deconvolved  from  the  regular  noise  associated 
with  an  imperfect  system.  The  noise  floor  of  the  instrument  was  too  high  for  the  percent 
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change  capable  of  the  lamp  combined  with  the  signal  magnitude  to  be  able  to  distinguish 


significant  SC  As  in  the  imaginary  spectrum.  Using  the  same  analysis  found  in  the  model 
is  difficult  when  the  overall  structure  of  the  imaginary  spectrum  is  hidden  within  the 
instrument  noise,  however,  over  certain  regions  in  the  speetrum  with  sharp  changes, 
imaginary  peaks  above  the  noise  level  eould  be  deteeted. 

Figure  17  provides  the  time  evolution  of  the  interferogram  ZPD  voltage  for  each 
scene  change  interferogram.  Eaeh  point  represents  a  different  interferogram  ZPD 
voltage.  Assuming  no  dwell  time  between  scans  and  using  a  scan  rate  of  1.17  Hz,  the 
starting  voltage  level  for  the  scene  change  interferogram  ZPD  was  assumed  to  be  the 
voltage  level  linearly  interpolated  to  3/7  seconds  before  the  ZPD  voltage. 
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Figure  17.  SCA  interferogram  ZPD  voltage  evolving  with  time. 
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Similarly,  the  ending  scene  change  interferogram  ZPD  voltage  was  determined  by 
linearly  interpolating  3/7  seconds  after  the  scene  change  interferogram  ZPD  voltage  on 
Figure  17. 

Figure  18  shows  the  spectrum  of  Measurement  10  at  17.14  seconds  which 
corresponds  to  the  beginning  of  the  decline  in  input  power.  The  spectral  features  from 
approximately  2500  cm  *  to  5500  cm  *  can  be  seen  to  represent  the  spectral  features 
caused  by  the  filament  in  the  lamp  but  they  were  unable  to  change  fast  enough  to  produce 
sizeable  SCA  magnitudes.  The  imaginary  components  rising  above  the  noise  level  near 
3900  cm  *,  4900  cm  *,  and  5250  cm  *  are  noticeable  in  the  spectra  taken  close  to  the  initial 
starting  point  of  the  radiance  decline.  The  3900  cm  *  and  5250  cm  *  peaks  could  be 
attributed  to  the  change  in  radiance  over  the  acquisition  time  being  able  to  create 
imaginary  spectral  features  that  could  rise  above  the  noise  floor  however  the  4900  cm  * 
peak  does  not  coincide  with  a  highly  structured  spectrum  so  its  cause  is  unknown.  The 
spectral  features  found  around  2000  cm  *  served  as  a  reference  for  every  spectrum 
because  they  did  not  significantly  change  during  the  entire  data  acquisition. 
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Figure  18.  Measurement  10  spectrum  at  17.14  seconds  with  a  ZPD  voltage  =  1.481 
V,  imaginary  on  top  and  real  on  bottom,  of  a  decaying  incandescent  lamp. 

Measurement  time  corresponds  to  a  few  SCA  spectral  peaks  in  the  imaginary 

spectrum. 

As  the  lamp  radiance  begin  to  decline,  spectral  features  were  no  longer  easily 
observable  because  the  overall  signal  magnitude  was  too  small.  Figure  19  shows  the 
spectrum  of  Measurement  25  which  corresponds  to  near  the  middle  of  the  210  W  to  0  W 
range  at  42.86  seconds.  The  radiance  magnitude  was  unable  to  provide  significant  and 
easily  recognizable  imaginary  spectral  features  as  seen  before. 
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Figure  19.  Measurement  25  spectrum  at  42.86  seconds  with  a  ZPD  voltage  =  0.47  V, 
imaginary  on  top  and  real  on  bottom,  of  a  decaying  incandescent  lamp. 
Measurement  time  corresponds  to  a  radiance  magnitude  and  change  that  created  no 
noticeable  SCA  imaginary  spectral  features. 

To  compare  the  magnitude  of  the  few  SCA  spectral  lines  in  the  scene  change 
spectrum  to  the  change  in  ZPD  voltage  that  occurred  during  the  scan,  the  difference  in 
ZPD  voltage  was  calculated  using  interpolated  values  based  on  Figure  17.  To  quantify 
the  magnitude  of  SCAs,  an  RMS  NESR  was  calculated  for  wavenumbers  next  to  the  SCA 
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location  where  there  was  only  noise.  The  RMS  NESR  was  calculated  by  selecting  an 
interval  near  the  region  of  interest  and  extracting  the  largest  positive  and  negative  peaks 
to  determine  a  peak  to  peak  NESR.  The  RMS  NESR  is  then 

RMSNESR  = 

6 

Eor  Eigure  18,  the  few  imaginary  peaks  above  the  noise  floor  near  3900  cm'^ 
reached  a  value  of  54.8%  above  the  RMS  NESR  for  ZPD  voltage  changes  of  2.52%. 
Since  Im{L^}  >  Re{L^}  it  is  not  unreasonable  to  assume  that  any  corruption  in 

the  real  spectrum  is  small  and  the  accuracy  in  the  spectrum  is  only  limited  by  the  noise 
level  of  the  instrument.  While  the  ZPD  voltage  change  was  higher  for  Measurement  25 
at  7.14%,  the  signal  magnitude  was  smaller  than  the  noise  so  no  peaks  emerged  above 
NESR. 


(36) 
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V.  Conclusions 


Interpretation  of  scene  change  effects  in  observing  detonation  flreballs 

The  nature  of  scene  change  on  FTS  provides  an  extra  layer  of  complexity  in 
interpreting  spectra  from  transient  events.  Using  an  approximation  given  by  Equation 
(21)  can  mathematically  describe  this  complexity  when  the  temperature  change  creates  a 
change  in  the  scene  radiance  described  by  Kick  et  al. 

Regardless  of  mathematical  representation,  quantifying  the  structure  and 
magnitude  of  SC  As  and  analyzing  Im{L^}  acquired  by  a  FTS  as  a  percentage  of  Re{L^} 

provides  insight  into  the  approximation  of  Re{L^}  asL^p^  and  the  resulting  level  of 

corruption.  This  is  useful  in  understanding  the  accuracy  of  using  FTS  on  detonation 
fireballs.  When  is  highly  structured,  the  magnitude  of  SCAs  increases  because  of  the 

increased  derivative  found  where  spectral  features  are  present.  The  ratio  of  SCA 
magnitude  to  the  total  spectrum  also  increases  as  temperature  changes  and  in  general, 
Im{4}  >  Re{L^}  -  L^pjj .  These  trends  were  independent  of  initial  temperature,  spectral 

resolution,  or  .  For  a  MOPD  =  0.5cm,  a  drop  in  temperature  of  3.3%  during  data 
acquisition  created  a  RMS  value  of  A3/JW I  sr/  cm^  /  cm~^  and  interpolating  to  a 
temperature  change  of  3.95%  gave  the  desired  5/JW  /  srlcrn^  I  cm~^ .  For  the  same 
conditions,  a  temperature  decline  of  26.7%  gave  the  same  RMS  radiance  value  of 
5//VT  /  sr  !  cn^  I  cm~^ . 

Observing  a  lamp  with  declining  input  power  provided  only  a  few  spectral 
features  in  the  imaginary  spectrum  that  were  larger  than  the  noise  floor  associated  with 
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the  instrument.  If  SCAs  are  below  the  NESR  for  the  detection  instrument  then  the 
accuracy  of  the  spectrum  is  only  limited  by  the  NESR  of  the  instrument.  The  few 
spectral  features  were  54.8%  above  RMS  NESR  for  a  ZPD  voltage  change  of  2.52% 
during  data  acquisition  that  lasted  6/7  seconds.  Since  the  imaginary  SCAs  were  hardly 
noticeable,  it  is  not  unreasonable  to  assume  Re{L^}  -  is  negligible. 

Recommendations 

Future  work  could  analyze  specific  spectral  features  to  show  the  corruption  on 
spectral  features  of  interest  rather  than  a  RMS  value  for  the  whole  spectrum.  This  would 
be  especially  important  when  focusing  on  CO2  and  H2O  bands  commonly  found  in 
detonation  fireball  spectra.  The  addition  of  higher  order  terms  in  the  analysis  would  also 
more  accurately  describe  scene  change  effects  on  radiance  mathematically.  While  there 
might  be  significant  corruption  to  specific  spectral  features,  the  overall  spectrum  can  still 
be  fit  to  a  Planckian  for  blackbody  temperature  determination. 
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